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Do the environmental conditions affect the dust-induced fragmentation 
in low-metallicity clouds?: Effect of pre-ionization and 
far-ultraviolet/cosmic-ray fields 



(N 

o 

(N 



O 

u 

6 

c3 



> 

oo 

o 
q 

o 



X 



Kazuyuki Omukai 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
omukai@tap.scphys.kyoto-u.ac.jp 

(Received ; accepted ) 
Abstract 

We study effects of the fully ionized initial state, or pre-ionization, on the subsequent thermal evolution 
of low-metallicity clouds under various intensities of the external far-ultraviolet (FUV) and cosmic-ray(CR) 
fields. The pre-ionization significantly affects the thermal and dynamical evolution of metal-free clouds 
without FUV/CRs by way of efficient HD formation. On the other hand, the pre-ionization effect on 
the thermal evolution is limited in very low-density regime for more metal-enriched clouds ([Z/H] > —4) 
or those under modest FUV 10~ 3 ) or CR field (> 0.1 of the present-day Galactic disk levels). In 
any case, for > 10 8 cm~ 3 , neither the initial ionization state nor the irradiating FUV strength affect the 
thermal evolution. The dust cooling is an important mechanism for making sub-solar mass fragments in 
low-metallicity gas. Since this fragmentation occurs at the temperature minimum by the dust cooling at 
> 10 1 
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ll0 cm 3 , this process is not vulnerable either to initial ionization state or external radiation. 
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Introduction 



The first generation of stars, or the so-called popula- 
tion III.l stars (pop III.l stars; O'Shea et al. 2008), are 
formed out of the primordial pristine gas with a small ion- 
ization degree (~ 10~ 4 ) left from the cosmic recombination 
(Peebles 1968; Sasaki & Takahara 1993). Within miniha- 
los with > 10 6 M Q , the collapse of primordial clouds is in- 
duced by the H2 cooling, which sets a characteristic mass- 
scale of dense cores at several 100M Q (Bromm, Coppi, <fc 
Larson 1999; Abel, Bryan, & Norman 2002). Inside these 
cores, protostars eventually form as a result of further 
graviational collapse by the H2 cooling (Omukai & Nishi 
1998; Yoshida, Omukai, & Hernquist 2008). Although the 
exact value of final stellar mass is still elusive, it is specu- 
lated that the stars grow massive (> several 10M©) owing 
to combined effects of high accretion rate and low opacity 
in primordial gas (Omukai & Palla 2003; Bromm & Loeb 
2004; McKee & Tan 2008; Hosokawa et al. 2011). 

Despite with the same primordial composition, stars 
formed under influences of radiative/mechanical feedbacks 
from pre-existing stars arc expected to have a different 
characteristic mass-scale from that of the pop III.l stars, 
and thus the terminology of pop III. 2 stars have been pro- 
posed to specify them (O'Shea et al. 2008). Among var- 
ious feedback mechanisms, e.g., photodissociation by the 
stellar far-ultraviolet radiation (Omukai 2001; O'Shea k, 
Norman 2008), the sweeping by supernova blast waves 
(Mackey, Bromm, & Hernquist 2003; Salvaterra, Ferrara, 
& Schneider 2004; Machida et al. 2005), one of the most 
well-studied is the case of star formation in relic HII re- 
gions of the first stars (O'Shea et al. 2005; Nagakura 



& Omukai 2005; Machida, Omukai & Matsumoto 2009). 
The massive first stars ionize the surrounding material 
by UV radiation (e.g., Kitayama et al. 2004). Except 
two ranges (10-40A/ o and 140-260M Q ) in mass, they 
end their lives rather quietly without supernova explo- 
sions (Heger & Woosley 2002; Umeda & Nomoto 2002). 
In those relic HII regions, the recombination proceeds in 
timescale ~ 10 s yrs and the gas condenses again to form 
next-generation stars (Yoshida et al. 2007). Such pre- 
ionized primordial-gas follows a different thermal evolu- 
tion from the pristine primordial gas of pop III.l star for- 
mation (Uehara & Inutsuka 2002; Nakamura & Umcmura 
2002). In the pre- ionized gas, more H2 (with molecular 
fraction ~ 10~ 2 ) than in the pristine gas (~ 10 -3 ) forms 
via the electron-catalyzed reaction, 



H+e^H +7, 
H+H- -> H 2 + e 



(1) 
(2) 

and thus its cooling lowers the temerature below the min- 
imum attainable in the pristine gas (~ 200K). In such 
a cold gas, HD fractionation is strongly favored via the 
exothermic (with energy difference AE/Icb = 464K) reac- 
tion 

H 2 + D+^HD + H+ (3) 

and the resultant HD cooling further lowers the temper- 
ature to a few tens K, close to the CMB temperature 
in the high-redshift universe, Tqmb = 27.3K(1 + z)/10. 
Those HD-cooling clouds fragment into dense cores of a 
few IOMq, which bear stars at their centers incorporating 
a large part of the material in the cores (Yoshida, Omukai 
& Hernquist 2007; McGreer & Bryan 2008). Namely, the 
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history of ionization in a low-density medium, or "pre- 
ionization" , has a large effect on the subsequent thermal 
evolution of the metal- free gas, and possibly changes the 
characteristic stellar mass. 

For metal-enriched clouds, the pre-ionization effect has 
not been explored comprehensively so far. Some previ- 
ous studies, though, indicate that it has indeed effects on 
the thermal evolution at least in some range of metallic- 
ity For example, in the initially un-ionized cases with 
metallicity [Z/H] > —3.5 1 , cooling by fine-structure tran- 
sitions of carbon and oxygen lowers the temperature be- 
low the minimum value attainable solely by the H2 cool- 
ing in the pristine gas (~ 200K at ~ 10 4 cm -3 ) (Omukai 
2000; Bromm & Loeb 2003). On the other hand, from 
the pre-ionized initial condition, the minimum tempera- 
ture attainable does not change so much from that in the 
primordial gas even for metallicity as much as [Z/H]~ —2 
(Jappsen et al. 2007). 

The pre-ionization effect in the metal-free gas is largely 
due to the HD formation and cooling. HD is also known 
to be quite susceptible to irradiation of FUV as well as 
ionization by the cosmic rays. For example, even in the 
pre-ionized gas, existence of only a weak FUV field pre- 
vents HD formation and its cooling (Yoshida, Omukai, & 
Hernquist 2007b; Wolcott-Green & Haiman 2010). On the 
contrary, a certain amount of cosmic-ray (CR) irradiation 
(with a few % of the energy density in the present-day 
Galactic disk) enables HD formation/cooling in an ini- 
tially un-ionized cloud (Stacy & Bromm 2007). 

In this paper, we study whether the pre-ionization af- 
fects the subsequent thermal evolution of low-metallicity 
clouds under a variety of intensities of the external FUV 
and CR fields. We also aim to clarify the density range 
where the thermal evolution is influenced by the exter- 
nal FUV and CR irradiation. The current prefered the- 
ory expects that transition from the Pop III massive star 
formation mode to the Pop II/I low-mass star mode is 
caused by the modification of thermal evolution due to 
the accumulation of metals in the star-forming gas, albeit 
the exact value of the critical metallicity is still in dispute 
(Schneider et al. 2002; Bromm & Loeb 2003; Dopcke 
et al. 2011). However, if the influences of environmen- 
tal parameters, e.g., pre-ionization or external radiation 
fields, remain until high density where the fragmentation 
and low-mass core formation occur, we also need to care 
about those parameters in studying fragmentation proper- 
ties of low-metallicity clouds and thus the problem would 
become highly complex. 

Our calculations here, fortunately, revealed that the 
pre-ionization alters the thermal evolution only in low 
density except cases with very low-metallicity ([Z/H] < 
-4) and under the low-FUV (< 10- 3 )/CR fields (< 0.1 of 
the present-day Galactic disk levels). In addition, with a 
given metellicity, all the temperature evolution tracks with 
diffrent FUV/CR intensities converge before the density 
range where efficient cooling by dust triggers the cloud 



fragmentation. This demonstrates that in studying evo- 
lution of the dense cores, we need not care about the past 
ionization history except in the limited circumstances. 

This paper is organized as follows. In Sec. 2, we de- 
scribe the method of calculation and input physics in our 
model. In Sec. 3, the results for temperature evolu- 
tion of low-metallicity clouds are presented and analysed. 
Finally, in Sec. 4, we conclude the paper by summarizing 
our findings. 

2. Method of Calculation 

We calculate the thermal evolution of low-metallicity 
clouds by using a one-zone model, which is based on that 
by Omukai (2000) and Omukai et al. (2005), but with 
some modifications, including (i) chemical reaction net- 
work and its rate coefficients are updated and (ii) effects 
of FUV and CR irradiation, specifically, heating by pho- 
toelectric emission of dust grains and by CR ionization, 
and photo- and CR-induced chemical reactions, are added. 
The physical quantities, i.e., temperature, density, chemi- 
cal fraction, etc., calculated in this model are regarded as 
those at the center of the cloud. The quantities calculated 
by this model are known to reproduce well the central evo- 
lution by hydrodynamical calculations for low-metallicity 
clouds (Omukai, Hosokawa, & Yoshida 2010) as well as for 
primordial clouds under FUV irradiation (Shang, Bryan, 
& Haiman 2010). 

Without strong magnetic or turbulent supports, the dy- 
namical collapse of self-gravitating clouds is well described 
by the Larson-Penston type similarity solution (Larson 
1969; Penston 1968). During the collapse, the cloud devel- 
opes the core-envelope structure with homogeneous cen- 
tral region of nearly the Jeans length in size and the en- 
velope where the density decreases with radius as oc r~ 2 . 
By mimicking the above property of the Larson-Penston 
solution, we assume the density increases by 

dp = p 
dt t co i ' 

where p is the density of t he gas, and the collapse timescale 
t coi = ^/37r/32(l — f)Gp where / is the ratio of pressure 
gradient to gravity and given as a function of the effective 
equation of state (see Omukai et al. 2005). The tempera- 
ture evolution is followed by solving the energy equation: 



(4) 



de 



dt P dt \p 



A r 



where the specific thermal energy is written as 
1 kT 



(5) 



(6) 



7ad - 1 prrin 

with the ratio of specific heat 7 a( j, the temperature T, and 
the mean molecular weight /1, and, on the right-hand side, 
p is the pressure 

pkT 



P 



1 [Z/H] = \og(Z/ZQ), where Z is the mass fraction of metals in 
the medium. 



(7) 



and A not is the net cooling rate per unit mass. The cooling 
processes included in A nGt are the radiative cooling by the 
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H Lya transition, fine-structure line transitions of [CII], 
[CI], and [01], and rotational (and also vibrational for H2) 
transitions of molecules H2, HD, CO, OH, and H2O. In 
calculating the H2 cooling rate, we also include the H2-e 
and H 2 -H + collisional transitions in addition to the H 2 -H 
and H2-H2 collisions since, in a highly ionized gas, col- 
lisions with electrons and protons can potentially domi- 
nate the H2 excitation (Glover & Abel 2008). The newly 
added transition rates are taken from Draine, Roberge, <fc 
Dalgarno (1983) for the H2-C collisions and Gerlich (1990) 
for the H 2 -H+ collisions, following Glover & Abel (2008). 
The cooling rates by the other species are calculated as in 
Omukai et al.(2010). In evaluating the photon trapping 
effects, we use the column density Ah = ^hAj, recalling 
that the size of the core is about a Jeans length Aj . The 
heating processes are those associated with the H2 forma- 
tion, photoelectric emission from the dust grains, and CR 
ionization. The heating rates by photoelectric emission of 
dust grains and by CR ionization are taken from Wolfirc 
et al. (1995). 

The chemical reactions among species composed of el- 
ements H, He, C, and O are solved. We consider the 
following 50 species as in Omukai et al.(2005): H, H2, e~, 
H+, H+, H+, H", He, He+, He ++ , HeH+, D, D+, D", 
HD, HD+, C, C 2 , CH, CH 2 , CH 3 , CH 4 , C+, C+, CH+, 
CH+, CH+, Gnj, CH+, O, 2 , OH, CO, H 2 0, HCO, 
2 H, C0 2 , H 2 CO, H 2 2 , 0+, 0+ OH+, CO+, H 2 0+, 
HCO+, 2 H+, H3O+, H 2 CO+, HCO+ and H 3 CO+. The 
included reaction network among those species are up- 
dated from Omukai et al.(2005), whose reaction net- 
work is largely based on the UMIST database of Millar 
et al.(1997) in addition to the primordial-gas chemistry 
from several sources (e.g., Abel et al. 1997; Galli k, 
Palla 1998). Our new network is based on the updated 
UMIST database of Woodall et al.(2007), supplemented 
by primordial-gas chemistry of Glover & Abel (2008). 
Specifically, we construct our chemical network as fol- 
lows. First, we consider all the primordial- gas reactions 
in Glover & Abel (2008) except those involving doubly 
dcutcrated species D 2 and D^~. For uncertain three-body 
H2 formation rate, we choose the value recommended by 
Glover (2008). Then, from a large number of the reactions 
listed in Woodall et al. (2007), we take those involving 
the above 50 species only, but excluding the collider re- 
actions, which become important only in higher density 
than our calculated range. For overlapping reactions be- 
tween the two compilations, we use the rate coefficient 
presented in Glover & Abel (2008). Primordial-gas re- 
actions not in Glover & Abel (2008) but in Woodall et 
al. (2007) are also considered. Finally, photodissociation 
of H 2 and HD is added (see below). Notable addition 
to our previous chemical model includes the photo- and 
cosmic-ray induced reactions as well as cosmic-ray induced 
photon reactions. The rate coefficients for CR-/photo- 
induced reactions by Woodall et al. (2007) are presented 
for those for the intensities in the Galactic disk, i.e., the 
CR primary ionization rate Cdisk = 1-3 x 10~ 17 s _1 , and the 
FUV intensity with respect to the Draine field Go = 1.7 
(Draine 1978), where the Habing parameter G is defined 



by the FUV flux in the range 6-13.6 eV normalized by 
1.6 x 10~ 3 crg cm -2 s™ 1 (Habing 1968). We thus rescale 
the coefficients in Woodall et al. (2007) for CR-induced 
reactions and photo-ionization/dissociation by multiply- 
ing the factor C/Cdisk and Go/1.7, respectively. This scal- 
ing corresponds to the assumption that the FUV spec- 
trum is the same as the Galactic field, which can be in- 
validated, for example, for the radiation field dominated 
by massive pop III stars with high (~ 10 5 K) effective 
temperature. The photo-chemical reaction coefficients de- 
crease exponentially with the visual extinction Ay owing 
to the shielding by dust grains. The visual extinction is 
related with the column density of the hydrogen nuclei as 
Ay = 5.3 x 10~ 22 -/VHCm -2 (Z/Z Q ) under the assumption 
of grains with the same optical properties as in the solar 
neighborhood (Bohline et al. 1978; Rachford et al. 2009). 
For the photodissociation of H2 and HD, and CO, in ad- 
dition to the dust shielding, the shielding by H2 becomes 
important. For H2 photodissociation, we use the reaction 
rate by Draine & Bertoldi (1996) 



fc H2P d =4.50 x 10- n / sh exp(-2.5Av)G s" 



(8) 



where the shielding factor / s h are taken from Wolcott- 
Grccn & Haiman (2010). For the HD photodissociation, 
we use the same expression as H2 (eq. 8), but the proper 
shielding factor for HD (Wolcott- Green & Haiman 2010). 
For the photodissociation of CO, we modify the reaction 
coefficients from Woodall et al. (2007) by including the 
shielding by H2. We use the fit to Lee et al. (1996)'s 
results by Hosokawa & Inutsuka (2006), which is within 20 
% accuracy for Ah 2 < 8 x 10 21 cm -2 although substantially 
worse at higher Ah 2 : 



fccopd = 2.0 x 10~ 10 cxp ( -3.5 A v - 



1.6 x 10 21 cm 



(G /1.7) a" 



We do not take the shielding of CRs into account. 
Although the CR shielding becomes important for column 
density > 29g cm -2 (Umebayashi & Nakano 1986), the 
column density does not reach such high and the shield- 
ing does not have a large effect in our calculations. 

The initial number density is taken at 0.1cm -3 , which 
is a typical value in the relic HII regions of the first stars 
(Yoshida et al. 2007). We consider two types of the initial 
conditions: (i) un-ionized condition: the initial tempera- 
ture Ti n i = 1000K, the ionization degree j/i n i(e) = 3 x 10~ 4 , 
molecular fraction yi n i(H 2 ) = 2 x 10 -5 , and (ii) pre-ionized 
condition : T ini = 8000K, y ini (e) = 0.5, 2/ in i(H 2 ) = 0. The 
un-ionized condition (i) is typical for the formation of the 
first stars (Yoshida et al. 2006), while the pre-ionized 
condition (ii) corresponds to the transitional state where 
the recombination is proceeding in the relic HII region 
(Yoshida et al. 2007a). The factor of a few variations in 
those initial chemical abundances do not alter the thermal 
evolution of gas presented below. 



4 



Omukai 



3. Results 

3.1. No FUV/CR Cases 

First, we see the cases with neither CR nor FUV field. 
Figure 1 presents the comparison of thermal evolution of 
low-metallicity gas in the un-ionized and pre-ionized cases. 

The lower panel shows the results in un-ionized cases. 
In low densities 10 4 cm -3 ), temperatures are almost the 
same for the clouds with [Z/H]< —4 as H2 is the domi- 
nant coolant in this density/metallicity range. For higher 
metallcity, i.e., [Z/H]> —3, cooling by the C and O fine- 
structure transitions breaks this degeneracy: the temper- 
ature becomes lower for higher metellicity. In the higher 
density range, temperature tracks deviate from the Z = 
one also for [Z/H] < —4, owing to cooling by H2 formed 
by the dust surface reaction, and by H 2 0, both of whose 
amounts are dependent on metallicity. 

Next we see the pre-ionized cases (upper panel). For 
[Z/H]< —4, HD dominates the cooling in < 150K, and 
thus the temperature tracks degenerate until ~ 10 8 cm -3 
where the gas is heated up to ~ 300K by the exothermic 
three-body H2 forming reaction. In comparison with the 
un-ionized clouds, which cool solely by H2, the tempera- 
ture in the pre-ionized clouds is significantly lower in the 
metallicity range [Z/H]< —4, thanks to the HD cooling. 
This effect, efficient HD cooling in the pre-ionized gas, is 
well-known for the primordial gas, i.e., the pop III. 2 star 
formation. For > 10 8 cm -3 , the warm environment makes 
H2 the more important coolant than HD, and thus the 
temperature tracks become the same as those in the un- 
ionized cases. With more metals ([Z/H] > —3), the cooling 
rate by the fine-structure lines exceeds that by HD. The 
temperatures are then similar in both un-ionized and pre- 
ionized cases, although slightly lower in the pre-ionized 
cases. In the case of [Z/H]= —1, the pre-ionized gas takes 
slightly longer time, i.e., higher density, to reach the CMB 
temperature of 30K. This just reflects the higher initial 
temperature in the pre-ionized case. 

To summarize, below the metallicity [Z/H]~ —3, HD is 
the dominant coolant for < 10 8 cm~ 3 in the pre-ionized 
gas, and the thermal evolution depends on the initial 
ionization state as is known for the pop III star forma- 
tion. For higher metallicity, on the other hand, the fine- 
structurc-linc cooling becomes more important and the 
pre-ionization has little impact on the subsequent ther- 
mal evolution. 

3.2. Cases only with FUV irradiation 

Next we see the cases with external FUV irradiation 
but without CR incidation. Panels in the Figure 2 corre- 
spond to the cases with different metallicitics Z = (left- 
top), [Z/H] = —5 (right-top), —4 (left-middle), —3 (right- 
middle), —2 (left-bottom) and —1 (right-bottom). In each 
panel, the solid and dashed lines indicate the un-ionized 
and pre-ionized cases, respectively, for five different FUV 
strengths, D = G /1.7 = (blue), 10~ 3 (green), 10~ 2 (or- 
ange), KU 1 (magenta) and I (red). In the pre-ionized 
cases (dashed) with strong FUV fields (Dq > 0.1), heat- 
ing by the photoelectric emission and cooling by the Lya 



emission balances, and the temperature remains almost 
constant at the initial value of 8000K for a while. In the 
un-ionized cases (solid), the temperature increases adia- 
batically from the initial 1000K to several thousand K, 
depending on the strength of FUV field, since the pho- 
todissociation prevents the H2 formation and its cooling. 
In both un-/pre- ionized cases, the temperature begins to 
fall down when the H2 self-shielding against the photodis- 
sociation becomes effective for [Z/H]< —3. With high 
enough metallicity [Z/H]> —2, the fine-structure-line cool- 
ing causes the temperature decrease before the H2 is self- 
shielded. 

As discussed in Sec. 3.1, without an FUV field, the 
thermal evolution of low-metallicity gas ([Z/H] < —4) bi- 
furcates according to its initial ionization state. This 
still holds true in weak enough FUV cases (see the cases 
with D < 1CT 3 ). With higher FUV intensity (D > 
10~ 2 ), however, the temperatures in un-ionized and pre- 
ionized cases quickly converge each other during the ini- 
tial temperature-decreasing phase. The dependence on 
the initial ionization state immediately disappears. This 
is because sufficient HD for cooling is not formed in either 
un-ionized nor pre-ionized case: even in the pre-ionized 
cases, H 2 photodissociation by FUV field prevents the gas 
to cool < 150K, thereby quenching HD formation, which 
proceeds only in the low-temperature environments. For 
[Z/H] ^ —3, as in no UV cases, the minimum tempera- 
tures are the same for pre- / un-ionized cases and no effect 
of pre-ionization is observed during the subsequent ther- 
mal evolution. 

It should be noted that, in all the cases, for > 10 8 cm~ 3 , 
neither the initial ionization state nor the irradiating FUV 
strength affect the thermal evolution. This is because the 
H2 fraction boosts up at ~ 10 8 cm~ 3 by the three-body re- 
action and its differences due to different initial ionization 
or FUV field strength in lower densities are wiped out. For 
high metallicity cases [Z/H] > —3, this convergence of the 
thermal evolution occurs at even lower density since the 
H2 cooling is only subordinate to the fine-strucutre-line 
cooling. 

3.3. Cases also with CRs 

Finally, we see the cases not only with FUV irradiation, 
but also with the CR incidation. In the solar neighbor- 
hood, the CR energy density is estimated, in terms of the 
primary ionization rate, to be (disk ~ 10 _17 s _1 (Webber 
& Yushak 1983), with about an order of magnitude un- 
certainty (e.g., Indrilo et al. 2007). Although the CR 
strength is, of course, far more uncertain in the early 
universe, about 1/100 of the Galactic value £ <~ 10 _19 s _1 
is claimed as an intergalactic CR field (Stacy & Bromm 
2007). In this section, we regard the CR intensity as a 
free paremeter and study cases with £ = 10 _19 ...10~ 17 s~ 1 . 

Figure 3 shows the weak CR cases of C = 10~ 19 s -1 . At 
first glance, these results resemble very well with those 
without CRs shown in Figure 2. However, at very low- 
metallicity [Z/H] < —4 and at the same time with low FUV 
intensity Dq < 10~ 3 , the temperature in the un-ionized 
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Fig. 1. Temperature evolution in clouds with six different metallicities (Z=0 (red), [Z/H]=-5 (magenta), -4 (orange), -3 (green), 
-2 (cyan), and -1 (blue)) against number density in the pre-ionized (upper panel) and for un-ionized (lower panel) cases. The cases 
with neither FUV nor CR incidation. The values of metallicities [Z/H] are indicated by numbers in the lower panel. The thin dashed 
lines indicate the constant Jeans mass, whose value is shown by numbers in the upper panel. 



cases is somewhat lower in the range of 10 5 — 10 7 cm 3 
than in the no CR cases, because the CR ionization en- 
ables the HD cooling (c.f. Stacy & Bromm 2007). The 
cases with higher CR intensity £ = 10 -18 and 10 _17 s _1 
are presented in Figures 4 and 5, respectively. The lat- 
ter CR value is as high as in the present-day Galactic 
disk. Most notably, with increasing CR strength, the 
pre-ionization effect becomes weaker. Without CR, the 
pre-ionization effect is clearly visible for clouds with low- 
metallicity [Z/H] < —4 and simultaneously with low FUV 
intensity Dq < 10 -3 (Sec. 3.2). With increasing CR inten- 
sity, the differences among the pre-ionized and un-ionized 
cases diminish even in those cases. Recall that, without 
CR, the HD formation/cooling is operative only in the 
pre-ionized environment. With some CR incidation, CR 
ionization induces the HD formation/cooling even in the 
un-ionzed condition. The differences in temperature evo- 
lution among pre-ionized and un-ionized cases thus be- 
come smaller for higher CR strength. Besides ionization, 
CRs also have the heating effect. This is clearly seen in 
Figures 4 and 5 as a rapid initial temperature rise in the 
un-ionzed cases with metallicity [Z/H]< —2. For density 
> 10cm -3 , this effect is not visible compensated by effi- 
cient cooling. 



4. Summary and Conclusion 

We have studied the effect of pre-ionization, as well as 
the far-ultraviolet irradiation and cosmic-ray incidation, 
on the thermal evolution of low-metallicity star-forming 
gas. Our findings can be summarized as follows: 

1. Without external FUV or CR irradiation, HD for- 
mation is promoted in the initially ionized ( "pre-ionized" ) 
clouds. This HD cooling makes the temperature in low- 
metallicity ([Z/H] < —3) clouds fall below the minimum 
value in the un-ionized gas attainable solely by H2 cool- 
ing (~ 200K). The temperature finally reaches a few 10 
K, similar to the high-redshift CMB temperature. With 
more metals, regardless of HD formation/cooling, the gas 
is able to reach such low temperature by the fine-structure 
line cooling. In this case, the pre-ionization does not have 
impact on the thermal evolution. 

2. Presence of a weak FUV field (G ~ 10 -2 ) pre- 
vents the HD formation/cooling even in the pre-ionized 
gas. This is because of the H2 photodissociation: with- 
out sufficient H 2 , the gas does not cool to < 150K, where 
the deuterium is rapidly converted to HD. In this case, 
the pre-ionization does not affect the subsequent thermal 
evolution of low-metallicity gas. 

3. The cosmic-ray ionization as high as £ ~ 10~ 18 s -1 in 
terms of the primary ionization rate, or ~ 10% of the solar 
neiborhood level, enables HD formation/cooling even in 
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Fig. 2. Thermal evolution of low-metallicity clouds with FUV irradiation, but without CR incidation. Panels correspond to the 
cases with different metallicitics Z = [left-top), [Z/H] = —5 (right-top), —4 (left-middle), —3 (right-middle), —2 (left-bottom) and 
— 1 (right- bottom). In each panel, the solid and dashed lines indicate the un-ionized and pre-ionized cases, respectively, for five 
different FUV strengths, Do = Go/1.7 = (blue), 10 — 3 (green), 10~ 2 (orange), 10 _1 (magenta) and 1 (red). The thin dashed lines 
indicate the constant Jeans mass, whose value is shown by numbers in the upper-left panel. 
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Fig. 3. Same as Figure 2 but with CR indication at the cosmic-ray primary ionization rate £ = 10 s 
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Fig. 4. Same as Figure 2 but for £ = 10 18 s 1 . 
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Fig. 5. Same as Figure 2 but for ( = 10 17 s 1 . 
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the un-ionized gas, thereby eliminating the pre-ionization 
effects even in very low-metallicity cases. 

In summary, the pre-ionization has remarlable effects 
only in clouds with very low mctallicity ([Z/H] ^ —3), 
and under low FUV (G < 10~ 2 ) or low CR (C < lO" 18 ) 
fields. 

In addition, 

4. the effects of external FUV/CRs on thermal evolu- 
tion vanishes until 10 8 cm~ 3 in all the cases we studied. 

The convergence of thermal evolution at high density 
has important consequence on the mass-scale of form- 
ing stars. In the framework of dynamical fragmentation, 
clouds fragments during efficient cooling phase, where 
the temperature decreases with density, while fragmen- 
tation hardly occurs in the temperature-increasing phase 
(Larson 1985, 2005; Li et al. 2003). The characteris- 
tic fragmentation mass-scales are thus set by the Jeans 
mass at the local minima of temperature. The low- 
metallicity clouds with — 5 ^ [Z/H] < — 1 has two tem- 
perature minima during the collapse; the first one in 
the lower-density (*~ 10 3 — 10 5 cm -3 ) due to the metal or 
H2, HD line cooling, and the second one in the higher- 
density (~ 10 10 — 10 15 cm -3 ) due to the cooling by dust 
thermal emission. Since the location of the first lower- 
density temperature minimum is vulnerable to the exter- 
nal FUV radiation, the fragmentation mass-scale owing 
to the line cooling can vary from cloud to cloud depend- 
ing on the intensity of surrounding FUV field. On the 
other hand, as the temperature evolution depends nei- 
ther on the pre-ionization nor external FUV radiation for 
> 10 8 cm -3 , the location of the temperature minimum by 
the dust cooling is robust. Thus, in studying the fragmen- 
tation by dust cooling we do not need to care about the 
external radiation. The results of numerical simulations 
for dust-induced fragmentation (e.g., Tsuribe & Omukai 
2006, 2008; Clark et al. 2008; Dopcke et al. 2011) remain 
valid also in the external radiation field. 

The author is very grateful for constructive comments 
by anonymous referee. This study is supported in part by 
the Grants-in-Aid by the Ministry of Education, Science 
and Culture of Japan (2168407, 21244021). 
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